# Read in both the oil spill and CA counties data
ca_counties <- read_sf(here("data","ca_counties","CA_Counties_TIGER2016.shp")) %>%
clean_names()
oil_spills <- read_sf(here("data","oil","ds394.shp")) %>%
clean_names()
# Change the CRS of the spills data to match the CRS of the counties data
oil_spills <- st_transform(oil_spills, 3857)
# Creating an interactive map
tmap_mode(mode = "view")
tm_shape(ca_counties) +
tm_fill("aland", palette = "BuGn") +
tm_shape(oil_spills) +
tm_dots()